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ABSTRACT 



Aims. We demonstrate the feasibility of determining parallaxes for nearby objects with the Wide Field Camera on the United Kingdom 
Infrared Telescope (UKIRT) using the UKIRT Infrared Deep Sky Survey as a first epoch. We determine physical parameters for 
ULAS J003402. 77-005206. 7, one of the coolest brown dwarfs currently known, using atmospheric and evolutionary models with the 
distance found here. 

Methods. Observations over the period 10/2005 to 07/2009 were pipeline processed at the Cambridge Astronomical Survey Unit and 
combined to produce a parallax and proper motion using standard procedures. 

Results. We determined n = 79.6 ±3.8 mas, fi„ = -20.0 ±3.7 mas/yr and yu^ = -363.8 ± 4.3 mas/yr for ULAS J003402.77-005206.7. 
Conclusions. We have made a direct parallax determination for one of the coolest objects outside of the solar system. The distance is 
consistent with a relatively young, 1-2 Gyr, low mass, 13 - 20 Mj, cool, 550-600K, brown dwarf. We present a measurement of the 
radial velocity that is consistent with an age between 0.5 and 4.0 Gyr. 

Key words. Astrometry - Stars: low-mass, brown dwarfs, fundamental parameters, distances 



1. Introduction 

The prototype T dwarf was discovered in 1995 as a companion 
to the nearby M dwarf star G1229 (Nakajima et al. 1995), and 
today we know of 155 T dwarfs (www.dwarfarchives.org as of 
10/2009). The majority of t hese are a result of th e near-infrared 
2 Micron All Sky Survey ( iSkrutskie et"ani2006l 2MASS) and 
the deep optical Sloan Digital Sky Survey dYork et al.l l2000i 
SDSS). Once discovered, significant efforts were undertaken to 
determine their d istances ( Vrba et al. 2004; Tinney et al. 2003; 
iDahn et aLll2002l) to map out the lower end of the Hertzsprung- 
Russell (H-R) diagram and to constrain models. However, in 
these surveys only a handful of objects have spectral types T7 
or later and more examples of intrinsically fainter and cooler ob- 
jects are crucial to understand both such objects and, in turn, 
exoplanets of similar masses. 

In 2005 t he UKIRT Infrared Deep Sky Survey 
dLawrence et all |2007, hereafter UKID SS) began with the 
Wide Field Camera dCasali et all l2007l WFCAM) large-field 
infrared camera. This survey is going about three magnitudes 
fainter than the 2MASS survey and has already revealed a 
number of extremely faint and cool T dwarfs (Lodieu et al.' 
2009; Burningham et al . 2008; Pinfield et al. 2008; Lodieu et al. 



20071). As part of the foUowup for the UKIDSS surveys we have 



started a program to determine parallaxes for the coolest objects 
being discovered. 
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The first and most obvious application of measured distances 
is to derive a luminosity which is used to populate the H-R di- 
agram and constrain models. Distances are also used to iden- 
tify unresolved companions via over-luminosities, and to pro- 
vide space motions which provide indications of an object's ori- 
gin and age. 

There is a large degree of uncertainty in current model anal- 
yses of the late-type T dwarfs, due to known inadequacies in 
the atmospheric models (e.g. Leggett et al. 2009). However, 
evolutionary models of cooling brown dwarfs are well under- 
stood, and show that all brown dwarfs older than about 200 Myr 
have a radius within about 20% of Jupiter's (e.g. B urrows et al.l 
2001). Hence there is a tight relationship between luminos- 
ity and effective temperature (Tes) for brown dwarfs (via the 
Stefan-Boltzman law). Parallax determinations combined with 
flux-calibrated observed spectra provide luminosities, and hence 
the best measurements of T^ff for brown dwarfs currently avail- 
able. Once Tgff is known, gravity and metallicity can be con- 
strained by comparing synthetic spectra to the observations (e.g. 
ISaumon et alj2007l) . lBurningham et"an (l2009l) use luminosity to 
determine T^ff for the T8.5 Wolf 940B, and show that an anal- 
ysis using synthetic spectra alone would have overestimated the 
temperature by ~100 K, for this 600 K object. 

The latest spectral type cuiTently defined is T9 and 
three such object s are c u rrently known: ULAS J003402.77- 
005206.7 (Warre n et alj l20 07bi herea fter ULAS 0034), 
CFBDS J0059 10.90-0 11401.3 (Delor meetalJ l2008l) and 
ULAS J133553.45 -hi 13005.2 (.Burningham et al.ll2008h . They 
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have effective temperatures in the range 500 to 700 K, and 
masses of 5 to 50 My corresponding to ages of 0.1 to 10 Gyr 
(e.g. Leggett et al. 2009). Here we present the first results of 
the UKIRT parallax program for ULAS 0034. In Sect. 2 we 
describe the observations and procedures for the program, in 
Sect. 3 we provide the results for ULAS 0034 and in Sect. 
4 we find the most consistent set of physical parameters for 
ULAS 0034 by combining available observations and models. 

2. Observations and Reduction Procedures 

2.1. Imaging 

The imaging observations were all made on the UKIRT 3.8 m 
telescope using the WFCAM imager This is the instrument be- 
ing used to carry out the UKIDSS surveys, and the calibrated 
images wer e taken directly from the pipeline for that survey 
jlrwin et al.l 12004) . All observations are carried out in queue 
override mode, allowing us to be very flexible in the scheduling, 
maximising the parallax factor and observing close to meridian 
passage. The observational sequence we adopt is very similar to 
the y-band imaging sequenc e in the Large Area Survey (LAS) 
component of UKIDSS (e.g. IWarren et al.ll2007" a'). A longer ex- 
posure time was used to provide an increased signal-to-noise for 
the targets, which are in general at the faint limit of the LAS. 
Specifically, we make observations in a 5 jitter (dithered) 3.2" 
cross pattern, and at each jitter position we make 4 exposures 
in 2x 2 micro-stepped positions of 1.5 pixels, where each expo- 
sure consists of 2 co-added 10 second images. The total expo- 
sure time is therefore 5x4x2x10 = 400 seconds. In average 
conditions this provides a signal-to-noise of 100 at 7 = 18. 

ULAS 0034 was discovered in a sequence of images taken on 
2005/10/04. We use the same pointing, filter and detector setup 
as the J image from that night in our parallax sequence. This 
allows us to use the discovery image as the first epoch in the 
parallax determination, and we reduce the time required for that 
determination by about one year (the time between the discov- 
ery observation and the recognition that an object warrants the 
investment required to determine a parallax). Nominally, to dis- 
tinguish between parallactic and proper motion we aim to have 
observations that cover at least three years, hence the saving of 
one year is significant. 

Parallax observations require precise astrometry. 
Unfortunately WFCAM has a large field-of-view with sig- 
nificant radial distortion, and as all four detectors are offset from 
the optical axis they are subject to large astrometric distortions. 
To minimize the differential astrometric distortion between the 
discovery and subsequent observations, the original pointing 
orientation is always used unless the object is close to the edge 
of a detector, in which case we move the object towards the 
centre. We work under the assumption that the astrometric 
distortion does not change for the duration of the observational 
program. This assumption forms the basis of all small-field, 
high-precision, ground-based parallax programs, even in tele- 
scopes where the focal plane is considered astrometrically flat, 
and is required to be true f or WFCA M because of the known 
large astrometric distortions (llrwin et al. 2004). 

During our parallax campaign WFCAM, a forward 
Cassegrain focus instrument, was regularly unmounted for sig- 
nificant periods so that the Cassegrain instruments could be used. 
To check that unmounting and remounting does not lead to sig- 
nificant changes in the astrometric map, a comparison of the ob- 
servations of the same cluster field for the nights 2005/06/20 
and 2005/09/04 were made, which straddles a period where 
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Fig. 1. Object position difference between two observations of 
the same cluster made on 2005/06/20 and 2005/09/04 as a func- 
tion of J magnitude. Plus signs are individual objects, crosses 
are median values. 



WFCAM was unmounted and remounted. In Fig. [T] we plot the 
differences in local x coordinates, as derived from the Cambridge 
Astronomical Survey Unit (CASU) pipeline discussed below. 
The first night was adopted as the reference system and the sec- 
ond night transformed using a six-constant model (translation, 
scale and rotation) derived from common objects. Assuming the 
nights were of equal quality the per observation error will be ap- 
proximately equal to the combined error divided by V2. From 
this comparison we estimate that the per observation precision 
for objects to 7 = 18 is around 10 mas in both coordinates in- 
cluding any contribution from the unmounting and remounting 
of WFCAM. An examination of binned residuals as a function 
of detector position did not reveal any systematic patterns. 

The observations made with WFCAM are processed and re- 
duced using a dedicated pipeline being run by the CASU. Since 
the centroid of the objects is our fundamental observation, we 
have extensively tested the results of this pipeline and we now 
describe these tests. 

First we have taken the images as they come out of the CASU 
pipeline (i.e. dark-subtracted, flat-fielded, systematic noise re- 
moved, sky-corrected, and finally registered and stacked) and 
found centroids two ways: (1) fitting a two-dimensional gaus- 
sian to the object's point source function (PSF) and (2) fitting 
a one-dimensional gaussian to the object's marginal distribu- 
tions above the sky background. The two-dimensional gaussian 
fit is what is used in the Torino Observatory Parallax Program 
(Smart et al. 1999, hereafter TOPP) and the marginal distribu- 
tion fit is what is suggested by Stone (1989) in the presence of 
a high background. We use these centroids and those from the 
CASU pipeline to compare eight ULAS 0034 observations over 
the period 2005 to 2008, after a six-constant transform (transla- 
tion, scale and rotation) to a common system. In Fig.|2]we plot 
the root-mean-square (rms) of the x coordinate from the CASU 
pipeline of common stars. The median rms is 18 mas and us- 
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Fig. 2. Root-mean-square of the x coordinates for stars in the first 
8 images of the ULAS 0034 sequence which spans 2.6 years, as 
a function of J magnitude. The centroids were all derived using 
the CASU pipehne. 



ing the other centroiding methods the rms was identical. In the 
y coordinate the median rms was slightly larger, 20 mas, but the 
difference between the three centroiding methods was also neg- 
ligible. 

We have also tried r ecombining the o riginal micro-stepped 
images using a drizzle dFruchter & Hoolo ,2002) routine rather 
than the CASU standard routine. In the standard routine the 4 
micro-stepped images are interleaved, i.e. combined assuming 
a higher spatial resolution, to produce a combined image with 
pixel size 0.2" (half the 0.4" physical pixel size). This produces 
a spiky PSF as the seeing changes between the four images. The 
alternative approach is to combine the counts in pixels from the 
four images in an underlying higher resolution grid with appro- 
priate weighting (i.e. drizzling). This gives a much smoother 
PSF and retains the signal-to-noise, while losing some of the 
resolution as adjacent pixels are correlated. Using this process 
caused the ULAS 0034 median rms to increase slightly from 18 
to 20 mas. 

As a result of these tests we have decided to use positions 
coming directly out of the CASU pipeline, and work under the 
assumption that our median error is around 18 mas. We assume 
the larger error found in the ULAS 0034 sequence compared to 
the two-epoch cluster test is due to the fact that we are com- 
paring a wider range of observing conditions over a longer time 
span. We note that the centroiding errors provided by the CASU 
pipeline are good indications of the errors found in our test. 

Once the (x,y) coordinates have been determined, the par- 
allax and proper motion of ULA S 0034 are determined u sing 
the methods adopted in the TOFF (ISmart et al.ll2003ll2007 l). We 
limit the reference objects used to only those within 5 arcminutes 
of the target, as these are sufficient for a transformation, and lim- 
iting the area of the detector being modeled also limits possible 
differential astrometric distortion. We have made two significant 
changes to the procedures used in the TOFF: 



(1) As we are observing in the /-band where the atmospheric re- 
fraction is small we expect the differential reddening correction 

DCR), to be negligible. This is in agreement with Tinnev et al] 
120031) and hence we do not apply any correction. Once we have 

data on more targets we will estimate the DCR and quantify the 

effect this could have on the target parallaxes. 

(2) Since we have the SDSS colours in this region we use the 
method of Ivezic eLal. (2008) to determine photometric paral- 
laxes for the anonymous reference stars and from those we cal- 
culate the correction from relative to absolute parallax. This cor- 
rection is typically less than 2 milliarcseconds. 

In the field there are also significant numbers of galaxies 
which can be used to correct from a relative to an absolute sys- 
tem. For the proper motions we found the mean of the galax- 
ies was less than the error of that mean and hence, within their 
error, these proper motions are on an absolute system. For the 
parallaxes we would need to u se the parameters of the galaxies 
as constraints ( lEichhornlll997h . future work will investigate this 
possibility. 



2.2. Spectroscopy 

We have also obtained new near-infrared spectra of ULAS 
0034 as well as of the template T8 dwarf 2MASSJ0415 
(B urgasser et al.ll2002h . for comparison. By spectral type T8 the 
CH4 and H2O molecular bands that define the T sequence are 
practically saturated. Therefore, despite the low effective tem- 
perature of ULAS 0034, 55 - 600K (see next s ection), relative 
to 750K for 2MASS J0415 (ISaumon et al.ll2007t) . at low resolu- 
tion the near-infrared spectral differences are reasonably subtle. 
For this reason we obtained medium-resolution spectra in the Y 
and J bands to search for any spectral differences identifiable at 
the higher resolution. One possibilit y is the appearance of ab- 
sorption Unes due to NH3 (iLeggett et al.ll2007h . 

The observations were made with the ISAAC instrument on 
the ESO VLT in service mode, during several different nights 
within the period July to September 2008. We used a 1" slit and 
the medium-resolution (MR) grating, providing wavelength cov- 
erage of 0.999 - 1.046/imin the F band, and 1.199 - 1.258//min 
the J band. Total integration times for ULAS 0034 were 160 min 
in the Y band and 120 min in the J band, and for the much 
brighter source 2MASSJ0415 were 40 min in the Y band and 
also 40 min in the J band. We previously obtained a spectrum 
with ISAAC in the H band, of similar resolution, and the data 
reduction is described in Warren et al (2007b). We followed the 
same procedures for reducing the Y and J spectra here. 

When smoothed to the lower resolution of the original 
GNIRS spectrum (from Warren et al., 2007b) we find the new 
spectra show good agreement, confirming in particular the addi- 
tional absorption for ULAS 0034 compared to 2MASS J0415, in 
the blue wing of the J band. Nevertheless a search for individual 
absorption lines failed to reveal any new spectral features. 

The new /-band spectrum of ULAS 0034 is of substantially 
higher S/N than the F-band spectrum, and we use it here to mea- 
sure the heliocentric radial velocity. The average seeing, 0.85", 
was a little less than the slit width, and determines the wave- 
length resolution of the spectrum, 3.5 x 10 FWHM. The 
use of a slit width larger than the seeing is not ideal for the mea- 
surement of a radial velocity (which was not the primary purpose 
of the observations), since any mis-centering of the source in the 
slit will introduce a systematic error. We quantify this error in 
the next section, when discussing the radial velocity. 
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Fig. 3. Observations and solution for ULAS 0034. The point at 
(0.03, 1.0) is the discovery observation on 2005/10/04 which 
had slightly lower signal-to-noise than subsequent observations, 
hence the larger error The observation without error bars was 
made in twilight and in poor seeing, and is of low quality, and 
was rejected from the solution by three-cr clipping. 

3. Results 

In Table [T] we report the astrometric parameters derived for 
ULAS 0034. In Fig.|3]we plot the observations and solution de- 
rived. 



Table 1. Solution parameters for ULAS 0034 and corresponding 
absolute magnitudes 



a, 6 (J2000.S) 


00:34:02.70, -00:52:07.7 


Base epoch (yr) 


2008.1312 


Absolute parallax (mas) 


79.6 + 3.8 


yu„ (mas/yr) 


-20.0 ± 3.7 


Us (mas/yr) 


-363.8 ± 4.3 


Relative to absolute correction (mas) 


1.33 


Duration of observations, (yrs) 


3.81 


observations, reference stars 


17, 113 


Absolute Magnitude in J' 


17.65 + 0.11 


Absolute Magnitude in K' 


17.98 + 0.12 



' Absolute magnitudes based on apparent magnitudes published in 
Warren et al. (2007b) and including both parallax and apparent mag- 
nitude errors. 



We measured the heliocentric radial velocity of ULAS 0034 
relative to 2MASSJ0415 by cross-correlating the two J-band 
spectra, finding a value -43.6 ± 5.4 km/s, i.e. ULAS 0034 is 
blue-shifted relative to 2MASS J0415. The quoted uncertainty 
accounts for the S/N of the two spectra, but not for the system- 
atic error, due to mis-centering of one or other object in the slit. 
We estimate, conservatively, that the accuracy of centering the 
object in the slit was as poor as 0.2". Then the systematic ve- 



locity error, for a single measurement, would be ~10km/s. We 
computed this value by establishing the shift of the centroid of 
the light captured by the slit, given the measured seeing, and as- 
suming a Gaussian profile. This eiTor dominates the error budget 
for the radial velocity, and it would be possible to obtain a more 
accurate velocity with a different instrumental setup. Combining 
in quadrature the systematic error for two measurements with the 
random error, we obtain a total uncertainty of 15 km/s. Taking 
the mean of the (much more accurate) heliocen t ric ve locities for 
2MASS J0415 given in IZapatero Osorio et al] (120071) . we con- 
clude that ULAS 0034 has a radial velocity of -1-7 + 15 km/s. 
Combining this with the measur ed proper motio ns, parallax and 
the velocity of the Sun from Dehnen & BinnevI il998.) we find a 
space motion, in km/s, of U=23 + 3 (radially inwards), V= -8 ± 6 
(in the direction of Galactic rotation), and W=-9 +13 (vertically 
upwards). 

Finally, we compare apparent J magnitudes for our 3 year 
parallax sequence and find an internal standard deviation of 0.02 
magnitudes and a systematic variation of less than 0.01. Hence 
ULAS 0034 does not appear to be photometrically variable in 
the near-infrared. 

4. Discussion 

In the ULAS 0034 discovery paper. I Warren et alj (12007^. they 
use near-infrared spectroscopy and mid-infrared photometry, 
with model atmospheres, to derive parameters of 600 < Teff K< 
700 and 4.5 < log g < 5.1 for ULAS 0034, corresponding to 
a mass range of 15 to 36 My and an age range o f 0.5 to 8 Gyr. 
The i nferred photometric distance is 14-22 pc. lDelorme et al.l 
(l2008i) use near-infrared spectra alone, with models, to estimate 
re ff ~ 670 K aiid log g ^ 4.6 fo r ULAS 0034 consistent with 
the lWarren et alj ( l2007b! ) analvsis.'Leggett et all ( 12009 '. hereafter 
LCS09) used a model comparison to a more complete spectral 
energy distribution which included mid-infrared spectroscopy to 
constrain the physical properties of ULAS 0034 further. They 
find that 550 < T^e: K< 600, 4.0 <logg < 4.5 and 0.0 < [m/H] 
< -hO.3. These values implied a distance of 13 - 16 pc, a mass 
of 5 - 20 My and an age of 0. 1 - 2 Gyr. The parallax determined 
here allows us to further refine these results. 

The measured distance of 12.6 + 0.6, rules out the low grav- 
ity (log g = 4.0), very low mass (5 - 8 My) and very young age 
(0.1 - 0.2 Gyr) LCS09 fits to ULAS 0034. We can check that 
the distance derived here is consistent with the other solutions 
by deriving the bolometric luminosity and comparing this to the 
family of Teff , log g solutions calculated by structural models for 
a given luminosity (e.g. Saumon & Marley 2008, their Fig. 4). 

We derived the bolometric flux by interpolating between the 
flux calibrated 0.8 - 2.3 /im GNIRS and 8 - 15 /zm IRS spec- 
tra presented in LCS09 using T^ff = 550 K, [m/H] = H-0.3 and 
Teff = 600 K, [m/H] = 0.0 log g = 4.5 synthetic spectra. The 2.3 
- 8 jum region contributes 33 - 43% of the total flux. The con- 
tribution beyond 15 fim was determined assuming a Rayleigh- 
Jeans tail to the distribution and T^ff = 550 K or 600 K; this 
region contributes ~8% of the total flux. The model interpola- 
tion and the Rayleigh- Jeans extrapolation were scaled by adopt- 
ing the distance measured in this work, and the structural model 
radius of 0.11 R© as per LCS09. We find a luminosity of log 
L/Lq - -5.96 ± 0.06, including the 4.8% uncertainty in distance 
and an estimated total flux uncertainty of 13%. The flux uncer- 
tainty is dominated by the estimated 13% uncertainty due to the 
model interpolation section; other sources of error are the uncer- 
tainty in the flux calibration of the observed spectrum (3%) and 
the assumption of a Rayleigh- Jeans tail (2%). The luminosity is 
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consistent with the T^s = 550 - 600 K, log g = 4.5 LCS09 so- 
lution for ULAS 0034, implying an age of 1 - 2 Gyr and a mass 
of 13 - 20 My. Note that the measured parallax rules out the 
possibiHty that ULAS 0034 is a 550 - 600 K binary system as 
the distance would have to be larger than any of the LCS09 solu- 
tions. Finally we note the velocity and age are co nsistent with the 
youn ger populations (< 4 Gyr) in Figs. 5 & 6 of iHolmberg et al.l 
( 1200 9) but outs ide the very young object (< 0.5 Gyr) velocity 
box of Fig. 6 in lZuckerman & Songl(l2004 . 



Table 2. Summary of physical parameters of ULAS 0034 



Distance (pc) 


12.6 ± 0.6 


U,V,W (km/s) 


23±3, -8±6, -9 ±13 


Teff (K) 


550 - 600 


log.? 


4.5 


Mass (My ) 


13-20 


Age (Gyr) 


1-2 
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The current UKIDSS followup program has 1 1 targets which 
will result in equally precise parallaxes by the end of 2009. 
Following the success for ULAS 0034 we recently added 20 
other faint T dwarfs which will have enough observational cov- 
erage for a parallax determination by the end of 20 11 . 



5. Acknowledgments 

RS and HRAJ would like to acknowledge the support of Royal 
Society International Joint Project 2007/R3. This research has 
benefitted from the M, L, and T dwarf compendium housed 
at dwarfArchives.org and maintained by Chris Gelino, Davy 
Kirkpatrick, and Adam Burgasser We would like to express our 
thanks to Luca Rizzi for scheduling help and all the observers at 
UKIRT who have carried out the observations used here. The 
United Kingdom Infrared Telescope is operated by the Joint 
Astronomy Centre on behalf of the Science and Technology 
Facilities Council of the U.K. SKL's research is supported by 
the Gemini Observatory, which is operated by the Association 
of Universities for Research in Astronomy, Inc., on behalf of the 
international Gemini partnership of Argentina, Australia, Brazil, 
Canada, Chile, the United Kingdom, and the United States of 
America. The ESO ISAAC spectra described here were obtained 
under programme 381 .C-0036(A). 



References 

Burgasser, A. J., Kirkpatrick, J. D., Brown, M. E., et al. 2002, ApJ, 564, 421 
Burningham, B., Pinfield, D. J., Leggett, S. K., et al. 2008, MNRAS, 391, 320 
Burningham, B., Pinfield, D. J., Leggett, S. K., et al. 2009, MNRAS, 395, 1237 
Burrows, A., Hubbai'd, W. B., Lunine, J. I., & Liebert, J. 2001, Reviews of 

Modern Physics, 73, 719 
Casali, M., Adamson, A., Alves de Oliveira, C, et al. 2007, A&A, 467, 777 
Dahn, C. C, Harris, H. C, Vrba, R J., et al. 2002, AJ, 124, 1 170 
Dehnen, W. & Binney, J. J. 1998, MNRAS, 298, 387 
Delorme, P, Delfosse, X., Albert, L., et al. 2008, A&A, 482, 961 
Eichhorn, H. 1997, Astron. Astrophys., 327, 404 
Fruchter, A. S. & Hook, R. N. 2002, PASP 114, 144 
Holmberg, J., Nordstrom, B., & Andersen, J. 2009, A&A, 501, 941 
Irwin. M. J., Lewis, J., Hodgkin, S., et al. 2004, in Society of Photo-Optical 

Instrumentation Engineers (SPIE) Conference Series, ed. P. J. Quinn & 

A. Bridger, Vol. 5493, 411^22 
Ivezic. Z., Sesar, B., Juric, M., et al. 2008, ApJ, 684, 287 
Lawrence, A., Warren, S. J., Almaini, O., et al. 2007, MNRAS, 379, 1599 
Leggett, S. K., Gushing, M. C.. Saumon, D., et al. 2009, ApJ, 695, 1517 



